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Abstract 

The possibility of the slab or stripe phase separation (alternating 
ferromagnetic highly- conductive and insulating antiferromagnetic lay- 
ers) is proved for isotropic degenerate antiferromagnetic semiconduc- 
tors. This type of phase separation competes with the droplet phase 
separation ( ferromagnetic droplets in the antiferromagnetic host or 
vice versa). The interaction of electrons with optical phonons alone 
cannot cause phase-separated state with alternating highly-conductive 
and insulating regions but it stabilizes the magnetic phase separation. 
The magnetostriction deformation of the lattice in the phase-separated 
state is investigated. 

1 Introduction 

In magnetic semiconductors there are two main mechanisms of the charge 
carrier self-trapping: magnetic and polaronic. The former corresponds to 
trapping of the electron by a microregion of a changed magnetic phase cre- 
ated by this electron. This effect is most energetically favored in antifer- 
romagnetic semiconductors where the self-trapping occurs in ferromagnetic 
microregions 1 ' 2 . On the other hand, the magnetic semiconductors are partly 
polar crystals. For this reason polaronic effects should exist in them. The 
elecron-phonon interaction is not so strong in them as in the completely 
polar NaCl-type crystals, and only it alone cannot change the state of the 
charge carriers drastically. In particular, the small polarons are impossible 
in the magnetic semiconductors 3 . Otherwise, they would not be degener- 
ate at heavy dopings, and the charge carriers would not realize the indirect 
exchange in them. But, being added to the magnetic mechanism of the 
self-trapping, the lattice mechanism can considerably increase the stability 
of the self-trapped state 4 . 

In degenerate antiferromagnetic semiconductors the magnetic self- trap- 
ping is a cooperative phenomenon. It manifests itself as a phase separation 
into the highly-conductive ferromagnetic phase and insulating phase as a 



result of the fact that all the charge carriers concentrate inside the ferro- 
magnetic phase 5 ' 6 . Both these phases are charged oppositely. Due to the 
Coulomb interaction between them the phase intermixing takes place. 

Two main types of such intermixing were considered in an isotropic crys- 
tal: the droplet one 6 and the slab one 5 . In the former case at relatively 
low charge carrier densities, droplets of a metallic ferromagnetic phase arise 
inside an insulating antiferromagnetic host, with all the charge carriers con- 
centrated inside these droplets. At larger densities, the highly-conductive 
ferromagnetic phase occupies the larger part of the volume, and the antifer- 
romagnetic phase is represented by insulating droplets inside it. The slab 
structure corresponds to alternating highly-conductive ferromagnetic and 
insulating antiferromagnetic plane layers. As a particular case, it includes 
the stripe structure with the ferromagnetic phase consisting of monolayers 
which separate the antiferromagnetic slabs. As is well known, the stripe 
structrure attracts much attention in connection with the manganites, nick- 
elates, and high T c superconductors nowadays. 

In the pioneering paper [5] only very crude estimates of the energy of the 
slab phase-separated state were carried out. The electric properties of such 
a phase-separated state should be highly anisotropic: it should be highly- 
conductive along the layers and insulating across them for any charge carrier 
density in the range where the phase separation is realized. In the subsequent 
paper [6] much more accurate calculations were carried for the droplet model. 
This model made it possible to explain the transition from the insulating 
state to the highly-conductive state with increasing charge carrier density 
as a result of the percolation of the ferro-electronic liquid, i.e. of changing 
the topology of the ferromagnetic phase from multiply- connected to simply- 
connected. Such a transition was observed experimentally in EuSe and in 
other materials (see Refs. [3,7]). 

The problem of the phase separation becomes again highly actual due to 
enhanced interest to the manganites. The coexistence of the ferromagnetic 
and antiferromagnetic phases in them was established very long ago 8 . There 
are new experimental data which evidence coexistence of these types of the 
magnetic ordering (e.g., Refs. [9-16]). Though sometimes these data were 
interpreted in terms of the canted antiferromagnetic ordering, one should 
keep in mind that it is unstable against fluctuations of the charge carrier 
density 3 ' 17 " 21 . For this reason one can assume that these properties are 
related to the ferro- antiferromagnetic phase separation. This makes it nec- 
essary to reinvestigate theoretically the phase separation in the degenerate 
antiferromagnetic semiconductors more accurately. In particular, it is nec- 
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essary to find out whether the slab structure is really impossible in them. 

The present paper is devoted to this problem. A more general model than 
in Refs. [5,6] is used here which takes into account not only the magnetic 
self-trapping but also the polarization that. For this reason essentially new 
results are obtained not only for the slab structures but also for the droplet 
ones. It will be shown that for some systems, both polarizable and nonpo- 
larizable, the difference in the energies of the slab and droplet structures is 
so small that it is below the accuracy of the calculation. Hence, one cannot 
state with certainty that the slab structure is impossible in the isotropic an- 
tiferromagnetic crystals. The small difference in the energies assumes that 
a small force can transform one structure into the other even if the latter is 
normally unstable. For example, if normally the droplet structure is stable, 
a small uniaxial deformation can stabilize the slab structure. 

Independently of the type of the structure, the lattice polaronic effects 
can essentially influence the electron-magnetic phase separation making it 
considerably more stable. They increase the size of the ferromagnetic re- 
gions in both droplet and slab structures. On the other hand, they increase 
the density of electrons in the ferromagnetic phase. As a result, the total 
volume of the ferromagnetic phase and the total magnetization of the crystal 
become decreased. Hence, the electron-lattice interaction, suppressing the 
ferromagnetic phase, hinders the metal-insulator transition for the droplet 
geometry and the transition to the uniform ferromagnetic state for both 
geometries. 

Especially large is contribution of the polaronic effect in the case when 
the polaronic attraction between the charge carriers almost compensates the 
Coulomb interaction arising due to the mutual charging of both phases. If 
the former exceeds the latter, the phase separation of the types discussed 
here becomes impossible. The condition for this the inequality eo/2 < eoo. 
where eo and e^ are the static and high- frequency dielectric constants, 
respectively. But such a situation will not be considered here. There are 
antiferromagnetic semiconductors, in which eo/2 only slightly exceeds e^, 
e.g., EuSe with eo = 9.4 and e^ = 5.0 [22], so that the investigation of the 
polaronic effects close the stability boundary is a quite actual problem. 

Recently some authors expressed the point of view that the purely pola- 
ronic phase separation is possible in the manganites without magnetic phase 
separation (e. g., [23]). This seems hardly possible as such a state would 
correspond to formation of regions consistsing of the Pekar polarons. But, 
as is well known, the Pekar polarons are impossible for realistic parameter 
values. This qualitative consideration was confirmed numerically using the 
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procedure described below. The situation can be different when the small 
polarons are more energetically favored than the large polarons. Accordong 
to Ref. [24], the small polaron condensation into a string is possible. But, 
as was already mentioned, the small polarons do not exist in the magnetic 
semiconductors. 

It should be noted that other types of the electron-phonon interaction 
can be essential for certain systems. In Ref. [25] the idea was advanced that 
the interaction between the electrons and Jahn- Teller phonons can cause 
the phase separation in the manganites. But this problem will not discussed 
here. In Ref. [26] the magnetostriction at the magnetic single-electron self- 
trapping was investigated. In the present paper the modulation of the lattice 
due to the interaction of the slab structure with the acoustical phonons 
will be investigated. The problem of magnetostriction at the electron self- 
trapping became actual now due to its experimental discovery [27, 28]. 

2 The Hamiltonian of the system 

The Hamiltonian of the system consists of the Hamiltonian of the s — d— 
model H s d, the free phonon Hamiltonian H p and the Hamiltonian H sp of 
the s-electron-phonon interaction. 

H = H s d + H sp + H p + He (1) 

H sd = t ^2 a g,<7 a g+A,<7 ~ A ^2( sS s)aa' a *ga a ST' 

where the s-electron operators, corresponding to the con- 

duction electrons or holes at the atom g with the spin projection a, s is the 
s-spin operator, S g that of the d-spin of the atom g, A the vector connecting 
the first nearest neighbors. The crystalline structure is assumed to be simple 
cubic, the <i-spin magnitude being S. As usually, the inequality t S> IS 2 
must be met as the hopping integral t is of the first order of magnitude, and 
the d — d- exchange integral / of the second order in the small d- orbital 
overlapping. 

The free phonon Hamiltonian for essential phonons is given by the stan- 
dard expression 

ff p = 5>q & q 6 q (h=l) (3) 
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where 6* , 6 q are the operators of the phonons with the wave vector q. 
The electron - phonon Hamiltonian can be presented in the form 



H sp = - ^2[ C q ex P( i( lS)bqa* sa a sa + H.C.] (4) 

In the long-wave approximation, the interaction constant with the longitu- 
dinal optical phonons is given by the expression 
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where u) is the longitudinal optical phonon frequency, V = Na 3 is the volume 
of the crystal, a is the lattice constant. For the longitudinal acoustical 
phonons in the deformation potential approximation the electron-phonon 
constant is 



/ k 2 

C « = tE ^2^MN " = " (5i,) 

where E\ is the deformation potential constant, M the mass of the unit cell. 

The term He describes the Coulomb energy of interaction of the s- elec- 
trons with each other and with ionized impurities. The charge of the latter 
compensating the charge of the charge carriers is assumed to be distributed 
uniformly over the sample (the jellium model). 



3 The variational procedure 

The reason for the magnetoelectronic phase separation in an antiferromag- 
netic semiconductor is the fact that the charge carrier energy is lower for 
the ferromagnetic ordering than for the antiferromagnetic that. If the charge 
carrier density is insufficient to make the entire crystal ferromagnetic, the 
charge carriers can concentrate in a portion of the crystal and make it fer- 
romagnetic. 

. Obviously, the ferromagnetic and antiferromagnetic phases are charged 
oppositely. To diminish the Coulomb energy, these phases should intermix. 
Two geometries will be considered. The first of them assumes that a crystal 
becomes separated into alternating layers of a highly-conductive ferromag- 
netic phase and the insulating antiferromagnetic phase 5 . The second is the 
droplet geometry of the phase separated state investigated in Ref . [6] . Only 
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the charge carrier densities below the percolation threshold will be consid- 
ered here which corresponds to the ferromagnetic highly-conductive droplets 
inside the antiferromagnetic host. 

Like in Refs [5,6], a variational procedure will be used here to find the 
ground-state energy of the phase-separated system. One of the variational 
parameters is the ratio of the volume of the antiferromagnetic phase to that 
of the ferromagnetic phase X. For the slab geometry the second variational 
parameter is the thickness of the ferromagnetic layers L with the antiferro- 
magnetic layers of a thickness XL (in the lattice constants a). Both numbers 
L and XL are assumed to be integer. Then, if the lattice is simple cubic, 
the electrons are concentrated in the ferromagnetic layers with dimension- 
less coordinates f x = nL(l + X) + g x , (1 < g x < L) where their number per 
atom v(x) is equal to + X) (vq is the mean number of s-electrons per 
magnetic atom, n an integer). Outside these layers the number of the con- 
duction electrons is zero. For the droplet geometry the second variational 
parameter is the radius of the ferromagnetic droplet R. 

In the ferromagnetic phase all the charge carriers are assumed to be 
spin-polarized. For W = 12t <C AS this condition is met for all the charge 
carrier densities. For W = 12t » AS it is met when AS > /i where ft is 
their Fermi energy. 

Two approaches will be considered here. The first one, sketched in Ref. 
[5] and used in Ref [6], is valid for large L or R . In this case, in addition to 
the standard bulk kinetic energy Ep, the surface electron energy Eg must 
be introduced which stems from spatial quantization of the electron levels. 
Just competition of the surface energy and the Coulomb energy ensures a 
finite size of the ferromagnetic regions. The second approach assumes that 
a s-electron is confined to a single atomic plane, i.e., that its motion is two- 
dimensional {L = 1). The ground-state total energy of the phase-separated 
system Eps (per unit cell) can be represented in the form: 

Eps = E k - + E c + E dd + E pol (6) 

For large L or R, the kinetic energy of the system can be presented in the 
form 

E k = E B + E s (7) 

Here the bulk kinetic energy of the charge carriers Ep is given by the ex- 
pression valid for relatively small v: 

Ep = -6tv +K(l+X) 2 / 3 , K = 0.6^(^)^0, M^) = i(67r 2 ^) 2 / 3 (8) 
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. where vq is the mean s-electron number per unit cell. 

The term E$ is the surface energy of the electron gas. Like in Ref. [6], it 
is calculated in the Born-Oppenheimer approximation, i.e. by expansion in 
1/[vqL{\-\-X)] powers. The Dirichlet boundary conditions for the s-electron 
wave function are taken. For the slab structure 

2*8^(1 + 
and for the droplet structure 

For the case of L = 1 one obtains the following expression for the total 
kinetic energy 

E\ = —Atu + 27Tfo(l + X)t (10) 

The term Eq in Eq (6) representing the Coulomb energy of the phase- 
and charge-separated state is calculated using the jellium model: the total 
compensating charge with the density (- eu^) modeling the charge of ionized 
donors or acceptors is distributed uniformly over the crystal. For the slab 
structure one has: 

(no) 

For the droplet structure the calculation is carried using the Wigner 
spheres of a radius Rw surrounding the ferromagnetic droplets. The total 
charge of the Wigner sphere is zero so that Rw = + X) 1 / 3 . Hence, 

d 27re 2 ^ 2 [2X + 3-3(l + X) 2 / 3 ] 
Ec = 5^ (U6) 

The change in the d — ci-exchange energy due to the formation of the 
ferromagnetic regions is given by 

Em = D = -zIS 2 (12) 

where z = 6 is the coordination number for a simple cubic lattice, / < 0. 
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4 The polaronic energy 



The polaronic energy E po i entering Eq. (6) is found from the effective lattice 
Hamiltonian H po i = H p + < H sp > where the symbol < ... > denotes 
averaging over the ground state of the electronic-magnetic subsystem. This 
corresponds to the adiabatic approximation in the polaron theory which 
is justified even for a weak electron-phonon coupling since the nonuniform 
electron distribution arises not due to the phonons in a selfconsistent manner 
but due to the magnetization. This distribution is given in advance with 
respect to the phonons. 

For the slab structure on obtains from Eq. (4): 

< H sp >= -£C q i/(/x)[& q exp(iqf) +H.C], (13) 

f,q 

f (/*) = + X ) E t/- " < l + X)L]@[n(l + X)L + L- f x ], (14) 

n 

where Q(z) = 1 for z > and &(z) = for z < 0. 
One obtains from Eqs (13), (14): 

< H sp >= -iNyN^l+X) Y,[C(q x , 0, 0)6 fe)0) oP(?x)p(? a! )+H.C.], (15) 

exp(iq x L)-l 
P{(lx) = expfe) - 1 ■ (16) 

a _ exp[ig a GL(l + x)\ - 1 

expfeL(l + x)]-l ' 10 

where 

G = i(i^) (18 » 

is the number of pairs of the ferromagnetic and antiferromagnetic layers, 
L(l + X) the dimensionless thickness of a pair of these layers. As G — > oo, 
the quantity P(q x ) is nonzero only for g m = 27rm/[a(l + -X")L] where m is 
an integer. Hence, 



(19) 



y = L(l + M) + 1 [cog(7r(1 + X)L) _ 1]; Z = y_i (ig) 
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where S(n, m) = 1 for n = m and otherwise. 

Using Eqs (16), (17) and diagonalizing the averaged phonon Hamiltonian 
Hpoi (3), (13 - 15) by the well known procedure of the phonon operator shift 

b m ^b m + Pm , m = NvoC T ^\ (20) 

m =(^jl'°'°)' ""W. 

one obtains the following expression for the cooperative polaron energy per 
unit cell: 

E voi--^T JL._ Y — (21) 

.,9 1 — cos27rm/(l + X) 



1 -cos27rm/[(l + X)L] 

The term with m = is omitted from Eq (21) as it is compensated by 
interaction of phonons with the ionized impurity charge. 

For the droplets such an accurate calculation of E po i will not be carried 
out for a reason pointed out in the end of this Section. It is sufficient to 
evaluate this quantity again using the approach based on the Wigner spheres. 
Then in the adiabatic approximation (e.g., Ref. [3]) 



E} = J D\v)d\ (22) 



where D(r) is the electric induction. Calculating this quantity for each 
Wigner sphere and multiplying it by their number, one comes to the expres- 
sion for Ep ol . 

To determine using the same simplified treatment, one should note 
that one can introduce an analogue of the Wigner spheres for the slab struc- 
ture, too. Namely, one can introduce the Wigner layers with the total 
zero charge, consisting of the highly-conductive layer with the x coordi- 
nates between —L/2 and L/2 and two adjacent insulating half-layers be- 
tween L(l + X)/2 and L/2 and symmetrical one. For them in the adiabatic 
approximation (22) one arrives to an expression for E^ ol similar to that for 
E* ol . For them both 

E P = -E C ^ (23) 
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The question arises about the accuracy of Eq (23). One observes that for 
the slabs one obtains the same expression directly from Eqs (21), (11a) if the 
quantities L and X tend to oo and to 0, respectively. The same occurs for any 
L as X — > oo. Really, in this case the sum over m in Eq (21) becomes equal 
to 2X 2 L 2 Q(2) where ((n) is the Riemann zeta-function. Hence, Eq (23) and 
its analogue are valid for large-scale phase separation with relatively small 
increase in the density or for very large volumes of the antiferromagnetic 
phase independenty of the size of the ferromagnetic layers. As follows from 
Eqs (lla,b) and (23), the polaron instability (Ec + E p < 0) should occur 
for 

| < too (24) 

If these conditions are not met, Eqs. (21) and (23) lead to different 
results. For example, if one takes eo = 5 and = 2.5, then one obtains 
by numerical calculations from Eqs (21) and (11. a): E„/Ec = 1.0008 for 
L = 1 and k = XL = 1000 or for L = 1000 and k = 1 in agreement with 
Eq (23). But for L = k = 1 one obtains a value of 2.4317 for this ratio. 
This means that for such a structure the polaronic effects are considerably 
stronger than for thicker ferromagnetic or antiferromagnetic layers. Hence, 
in principle, they can stabilize the stripe structure when the ferromagnetic 
layers are reduced to monoatomic planes (the long-rang ferromagnetic or- 
dering is possible in the two-dimensional systems even at T / due to the 
magnetic anisotropy). 

As is seen from Eqs (23) and (lla,b), in the adiabatic approximation 
used here the polaronic energy vanishes for the uniform ferromagnetic state 
(for the slab structure it follows from the fact that according to Eq (21) 
p(q m ) = for X = 0). 

Now we return to the problem of a calculation of the polaronic energy 
in the manner similar to that leading to Eq. (21). As follows from above- 
said, such a calculation could be very useful for monoatomic sublattices 
which are the limiting case for the droplet model. But such sublattices 
have nothing in common with the phase separation as single atoms cannot 
form a ferromagnetic phase. Thus, this problem is beyond the framework 
of the present paper, and for this reason is not considered here. The Eq. 
(23) is sufficient for the droplet model which assumes that the size of the 
ferromagnetic droplets is large compared with the lattice constant. 
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5 Calculations and their results 



Now results of calculations will be presented. First, it should be noted that 
a necessary condition for a phase-separated state to be energetically favored 
consists in the requirement that the energy (6) be less than the energy of 
the uniform antiferromagnetic state E^p: 

E AF = -6tu + 2~ 2/3 K for W > AS (25) 

The latter Equation for the double-exchange system is obtained using an ex- 
pression for the charge carrier energy in a system with the double exchange, 
if the d-spins are quantum [3]. It should be noted that the term AS/ 2 does 
not disappear for any type of the magnetic ordering playing the part of an 
additive constant in the case of the double exchange. 

According to Eq (23), if one wishes to account for the polaronic effects, in 
the case of sufficiently large radii R or layer thicknesses L one can introduce 
the effective dielectric constant k, defined by the relationship 

I = ±_I = I_2. (27) 

k e e* e 

In this case one can carry out the minimization of the energies (6-9a,b ),(11 
a,b), (12) with respect to R or L in an explicit form, and then only the min- 
imization with respect to X should be carried numerically. If one interested 
only in the pure magnetic self-trapping, one should leave Eqs (lla,b) as they 
stand, and if one wishes to take into account the polaronic contribution, one 
should replace eo by k (27) in them. 

In this case one can exclude the quantities R ot L from the expressions 
(6) for the energy and make it a function of the only variational parameter X. 
For the droplets structure with both ferromagnetic and antiferromagnetic 
droplets this procedure was carried out in Ref. [6]. For the slab structure 
this program was not carried out yet. 

For the droplet model, as in Eq (6) only the terms Eg and Eq are in- 
dependent, one should replace the sum Eg + E c by its minimal value with 
respect to R: 

E d sc \.2v^[2X + 3 - 3(1 + X) 2/3 } 1/3 (1 + X) 2/9 (28) 
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e a 

and substitute the following expression for the radius: 

R= { 75^(1 + X)V3 6oa |V3 

\6 1 /3647r2/3 l/ 2 e 2[2X + 3-3(l + X)2/3] J { 



A similar calculation for the slab structure leads to the expressions 

2/3 / 2 \ 1 / 3 

2 V 3/ \ 16 



E , c - « (l)" (£f' 3 ^ + (30, 



L 



2 2 / 3 15(l + X)V3 eoa ] 1 / 3 



1 Q 7r 2/3 e 2 zy 2 X 2 



Further calculations should be carried numerically. 

Now some results of them will be presented. First, we consider the case 
when the the system is close to the boundary of the polaron instability. For 
this aim the same values of the dielectric constants will be taken as in EuSe 
[23]: eo = 9.4, = 5. As EuSe is an isotropic metamagnet with a very 
low critical field of about 2000 Oe, a value of D = 0.0001 is chosen (here 
and below all the typical energies are in the eV units). Other quantities 
are: t = 0.2, A = 1.2/7, e 2 /a = 2.7, S = 7/2 . Not all they coincide with 
the corresponding quantities in EuSe though are likely to be close to them. 
Nevertheless, the inequality W » AS is certainly met here, and the energy 
of the system will be counted off from the expression (25). 

The results of the calculations are depicted in Fig.l for both droplet and 
slab structures. Here and below the upper indices p and r will denote the 
polarizable and rigid lattices, the lower indices d and s the droplet and slab 
structures, respectively. The parameters X?, X^, BP and L p are presented 
in Fig.l For the droplet model the calculations are carried only down to 
the X p value of 1, which corresponds to the electronic and ferromagnetic 
percolation. Calculations of the polaronic effects for X p < 1 should be 
carried out using another procedure. 

As is seen from Fig. 1 , the parameters X v d and X p decrease and hence the 
ferromagnetic portion of the crystal increases with increasing charge carrier 
density uq. Respectively, the sizes of the ferromagnetic regions HP and LP 
increase. For the slab structure an abrupt concentration phase transition 
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into the uniform ferromagnetic state takes place though the X p jump is 
rather small. As the polaronic contribution to the energy is small at X p — > 0, 
this effect is not related to the polaronic effect. A corresponding calculation 
for the unpolarized droplet model using results of Refs [3,6] shows that the 
phase separation also disappears abruptly. 

Comparing numerical results for these two models, one should point out 
that their energies and X p values are close to each other. In both cases the 
energy per conduction electron decreases smoothly with increasing electron 
density. For the droplet structure this energy E d changes from - 0.2089 at 
i/ = 4 xlO- 4 (X p d = 5) to - 0.2245 at vq = 1.2 xl0~ 3 (X p d = 1). For the 
slab structure this energy E P is equal to - 0.2058 at u = 4xl0 -4 (X P = 
4.6), - 0.2256 at u = 1.2 xlO" 3 (X P = 1), and - 0.2387 at v = 1.98 xl0~ 3 
(Xf = 0). It is interesting to note that the densities corresponding to equal 
volumes of the antiferromagnetic and ferromagnetic phases coincide for both 
models. 

One should compare results obtained with accounting and without ac- 
counting for the lattice polarization. It is advisable to carry out this compar- 
ison for the density v$ = 0.00095 corresponding to the percolation threshold 
in the rigid lattice. One obtains: for the droplet model in the rigid lattice 
the energy per electron E r d = - 0.2054, X d = 1.0, R r d = 10, the number of 
electrons in the droplet N d = 8,3. The small parameter of theory ensuring 
applicability of Eq. (9a) s r = l/[(6vr 2 ) 1 / 3 iy 1/3 (l + x) 1 / 3 ] is 0.21 which justi- 
fies use of the many-electron approach adopted here. With accounting for 
the the polaronic effects one has: the energy per electron E d = - 0.2202, x p d 
= 1.6, R p = 16, N p d = 44, s p = 0.12. 

Obviously, the polarization increases the number of the electrons in each 
droplet by a factor of 4.75, whereas its volume is increased only by a factor 
of 4.1. Hence, the polarization increases the density of conduction electrons. 
The total magnetization of the crystal is proportional to the relative volume 
of the ferromagnetic phase 1/1 + X. Thus, the polarization decreases the 
total magnetization by about 30 %. Simultaneously, the system is shifted 
from the percolation threshold to the insulating side. 

Now these results will be compared with those obtained for the slab 
model with the same vq value of 0.00095. Without polarization, one obtains: 
E r s = -0.2073, X r s = 1.0, the ferromagnetic layer thickness U = 8.5. With 
polarization, E p = -0.2208, X P = 1.5, L p = 13. Though, formally, the energy 
of the slab structure is less than that of the droplet structure, the difference 
between then amounts only to 0.3 % which is certainly considerably less than 
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the accuracy of the present calculations. For this reason it is impossible to 
decide which structure should be realized in reality. The reduction of the 
magnetization caused by the lattice polarization in the slab structure is close 
to that in the droplet structure. 

Now we consider a double-exchange system of the manganite type with 
values of the dielectric constants coinciding with those of the La- based 
manganite: eo = 5,6^ = 3.4 [29]. The s — d- exchange integral A is assumed 
to be very large (its exact value is inessential as its enters the expression for 
the energy only as an additive constant ), and the following values will be 
taken: t = 0.2, D = 0.01, S = 2. Then one finds for the density u Q = 0.015 
close to the percolation density: for the unpolarized droplet structure the 
energy E d = - 0.06 is minimal for X d = 1.1 with R r = 4.2, N d = 10 and s r 
=0.19 (the energy is counted off from the energy U (26)). If one takes into 
account the polarization, then E p = -0.09, R p d = 4.8, N p d = 16, s p = 0.17 
with X p d = 1.3. 

In the case of the slab structure for the same values of parameters one 
finds: for a rigid lattice E r s = - 0.0708 with X r s = 1.1 and U = 3.44, for a 
polarizable lattice E p = -0.0986 with X P = 1.3 and L p = 3.84. 

As is seen from these data, the main regularities established above for 
the case of the wide bands W S> AS remain in force for the case of the 
double exchange. The difference in the energies of the slab and droplet 
structures is larger in this double-exchange system case which points to the 
droplet structure being preferential. One can also give examples where the 
situation is opposite. Being unable to predict in some cases which structure 
should be realized, I would like to point out that just due to this small 
energy difference a small external force can cause transition of one structure 
to the other. For example, if the droplet structure is normally realized then 
a small uniaxial stress increasing the hopping integral t in the corresponding 
direction can stabilize the slab structure. 

It is still necessary to discuss the possibility of the ferromagnetic stripe 
structure (L = 1). It is clear that the single layer ferromagnetic phase 
could be possible only in the double exchange systems since in the opposite 
case W 3> AS the energy of the two-dimensional motion of the electron 
is by 2t = W/6 higher than in the three-dimensional case. Meanwhile, 
the gain in its energy due to the ferromagnetic ordering is only AS/ 2. On 
the other hand, as was shown in the preceding section, the polaronic effect 
is enhanced in the structures with alternating ferro- and antiferromagnetic 
atomic planes. For this reason one can expect the existence of such structure 
in the double-exchange systems. 
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I managed to prove the possibility of such structure in the limit S —> oo. 
The following values of the parameters were used: eo = 5, eoo = 3.4, t = 
0.1, D = 0.01. Then one finds from Eqs (6) and (10) that for v Q = 0.02 
the system with L p = 1 is the most energetically favored for k = 3 (3 
antiferromagnetic layers per ferromagnetic layer). For v§ = 0.05 the most 
energetically favored structure is L p = 1, k = 1 with alternating ferro-and 
antiferromagnetic atomic planes. For small S, I did not find parameter 
values at which the stripe structure is most energeticlly favored, i 

It remains only to consider the magnetostiction deformation of a sample 
in the case of the slab phase separated state. For this aim one should use 
Eqs (20), (5b) yielding the shifts /3 q of the operators for the longitudinal 
acoustical phonons. Going over from the phonon operators to the normal 
coordinates in the standard manner, one obtains the following expression for 
their shift 5Q q : 

SQ « = {^M^ + ^ < 32 > 

where M is the mass of the unit cell. Using this Eq, one obtains for the 
displacement of the atom g along the x axis related to the wave with the 
vector m (20): 

u g = ,4 m exp(img), A m = -^=8Q m 

V N 

An _ Eim\p m -p* m \vp . . 

As follows from Eq (33), the lattice becomes modulated not only with the 
periodicity of the phase-separated structure but also the higher harmonics 
are present. 
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Captions 

Fig.l. The ratio X p of the volumes of the antiferromagnetic and fer- 
romagnetic phases for the droplet (X%) and slab (X p ) structures in the 
polarizable lattice, as well as the radius R p of the ferromagnetic droplet for 
the former and the thickness L p of the ferromagnetic layer for the latter, 
as functions of the relative average charge carrier density v$. R p and LP 
are expressed in the lattice constants. Numerical values of the parameters 
used, which correspond to W 3> AS, are presented in the text just below 
Eq. (31). 
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